Bistability and chaos at low-level of quanta 
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We study nonlinear phenomena of bistability and chaos at a level of few quanta. For this pur- 
pose we consider a single-mode dissipative oscillator with strong Kerr nonlinearity with respect to 
dissipation rate driven by a monochromatic force as well as by a train of Gaussian pulses. The quan- 
tum effects and decoherence in oscillatory mode are investigated on the framework of the purity of 
states and the Wigner functions calculated from the master equation. We demonstrate the quantum 
chaotic regime by means of a comparison between the contour plots of the Wigner functions and the 
strange attractors on the classical Poincare section. Considering bistability at low-limit of quanta, 
we analyze what is the minimal level of excitation numbers at which the bistable regime of the sys- 
tem is displayed? We also discuss the formation of oscillatory chaotic regime by varying oscillatory 
excitation numbers at ranges of few quanta. We demonstrate quantum-interference phenomena that 
are assisted hysteresis-cycle behavior and quantum chaos for the oscillator driven by the train of 
Gaussian pulses as well as we establish the border of classical-quantum correspondence for chaotic 
regimes in the case of strong nonlinearities. 

PACS numbers: 42.65.Pc, 42.50. Dv, 05.45.Mt 



I. INTRODUCTION 

Nonlinear dissipative oscillator (NDO) operated in the 
quantum regime is becoming significant in both funda- 
mental and applied sciences, particularly, for implemen- 
tation of basic quantum optical systems, in engineering 
of nonclassical states and quantum logic. The important 
implementations have recently been realized in the con- 
text of superconducting devices based on the nonlinearity 
of the Josephson junction (JJ) exhibiting a wide variety 
of quantum phenomena (see, the Reviews [1], |2] and [3]- 
[9]). In some of these devices, dynamics are analogous 
to those of a quantum particle in an oscillatory anhar- 
monic potential [10) . The single nonlinear oscillator and 
systems of nonlinear oscillators also consist of basic the- 
oretical models for various nano-electro-mechanical and 
nano-opto-mechanical devices. In the last decade there 
was exciting technological advances in the fabrication and 
control of such devices [ll]-[14] that are attracting inter- 
est in a broad variety of research areas and for many 
possible applications due to their remarkable combina- 
tion of properties: small mass, high operating frequency, 
large quality factor, and easily accessible nonlinearity. 
Nanomechanical oscillators are being developed for a host 
of nanotechnological applications. They are ideal candi- 
dates for probing quantum limits of mechanical motion 
in an experimental setting. Moreover, they are the ba- 
sis of various precision measurements |15j-|17j. as well as 
for basic research in the mesoscopic physics of phonons 
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[IH] , and the general study of the behavior of mechanical 
degrees of freedom at the interface between the quantum 
and the classical worlds [19]. 

The efficiency of quantum oscillatory effects requires a 
high nonlinearity with respect to dissipation. However, 
for weak damping even small nonlinearity can become im- 
portant. Development of driven NDO in quantum regime 
requires cooling these systems to their ground state and 
significant advances have been made in cooling the sys- 
tems to far below the temperature of the environment 

It is well known that classically, the nonlinearity- 
induced dependence of the oscillatory frequency or the 
amplitude usually leads to bistability of the driven NDO 
[22] ■ The corresponding systems describe amplifiers 
which are ubiquitous in experimental physics. The 
Josephson junction amplifier has been discussed in [2Sj 
as a bistable amplifier which has particular interest. The 
dynamical bifurcation of a rf-biased Josephson junction 
was proposed to be used as a basis for the amplification 
of quantum signals in [211 130] ■ Though the bistability 
has been usually treated as a classical signature of the 
NDO the quantum dynamics in the bistable region has 
been a new subject in the past years [5T|-|34j. 

Quantum dynamics of an oscillator is naturally de- 
scribed by Fock states, that have definite numbers of 
energy quanta. However, these states are hard to cre- 
ate in experiments because excitations of oscillatory sys- 
tems usually lead to the production of coherent states but 
not quantum Fock states. Nevertheless, quantum oscilla- 
tory states can be prepared and can be manipulated by 
coupling oscillators to atomic systems driven by classi- 
cal pulses. The systematic procedure has been proposed 
in Ref. [35] and has been demonstrated for determin- 
istic preparation of mechanical oscillatory Fock states 
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with trapped ions [3B], in cavity QEDs with Rydberg 
atoms [37] and in soUd-state circuit QED for determin- 
istic preparation of photon number states in a resonator 
by interposing a highly nonhnear Josephson phase qubit 
between a superconducting resonator |38| . 

For the NDO in quantum regime the nonhnearity 
makes frequencies of transitions between adjacent osciUa- 
tory energy levels different. Thus, strong nonhnearity en- 
ables spectroscopic identification and selective excitation 
of transitions between Fock states. Thus, it is possible in 
this regime to prepare the NDO at low- level of quanta. In 
this approach, it has been shown that the production of 
Fock states, Fock states superpositions or qubits can also 
be realized in over-transient regime of an anharmonic dis- 
sipative oscillator without any interactions with atomic 
and spin-1/2 systems and with complete consideration of 
decoherence effects [39]. For this goal the strong Kerr 
nonhnearity as well as the excitation of resolved lower 
oscillatory energy levels with a specific train of Gaussian 
pulses have been considered. 

In this paper we consider the NDO in the regime of 
low-level of excitation for the study of the problems of 
quantum bistability and chaos. Therefore, the goal of 
the paper is twofold. In one part, we consider the bista- 
bility on a few oscillatory excitation number considering 
the NDO driven by monochromatic force. The bistabil- 
ity on a few excitation number is attractive for ultra-low 
power operation, but it has practical problems related 
with quantum fluctuation-induced spontaneous switch- 
ing. In this part, we also demonstrate the production of 
quantum interference between bistable branches for the 
NDO driven by a train of Gaussian pulses. 

The other part of the paper is devoted to investiga- 
tion of quantum chaos in low-level excitation regime of 
the pulsed NDO. Much research on the subject of clas- 
sical and quantum chaos has been done on the base of 
the kicked rotor, which exhibits regions of regular and 
chaotic motion. Its experimental realization, and obser- 
vation of the models dissipation and decoherence effects 
are carried out on a gas of ultracold atoms in a magneto- 
optical trap subjected to a pulsed standing wave [40] . 
[JT]. In Ref. [12] it was proposed to realize the para- 
metrically kicked nonlinear oscillator model in a cavity 
involving Kerr nonhnearity. It was also shown that a 
more promising realization of this system, including the 
quantum regime, is achieved in the dynamics of cooled 
and trapped ions, interacting with a periodic sequence of 
both standing wave pulses and Gaussian laser pulses [13] . 
In addition to these and the other important input (see, 
also [44j . |45| ) in this paper we consider quantum dissi- 
pative chaos of the NDO driven by a train of Gaussian 
pulses in a strong quantum regime and in complete con- 
sideration of dissipation and decoherence. In this way, 
quantum dissipative chaos at the limit of low-level of ex- 
citation numbers will be considered. 

The paper is arranged as follows. In Sec. II, we shortly 
describe a pulsed NDO. In Sec. Ill, we study bistability 
at a level of few quanta for the NDO driven by monochro- 



matic force. In Sec. IV, we consider the production of 
quantum interference for bistable regime of the NDO 
driven by the train of Gaussian pulses. In Sec. V, we 
study quantum dissipative chaos at limit of low-level of 
excitation numbers for the NDO driven by the train of 
Gaussian pulses. We summarize our results in Sec. VI. 

II. THE MODEL: SHORT DESCRIPTION 

The Hamiltonian of an anharmonic-driven oscillator in 
the rotating-wave approximation reads as: 

H = hl\a)a + hx{a)af + hf{t){na) + Wa). (1) 

Note that time dependent coupling constant ^f{t) that is 
proportional to the amplitude of the driving field consists 
of Gaussian pulses with duration T which are separated 
by time intervals r as follow 

/(t) = ^e-(*-*o--)VT^. (2) 

Here, a\ a are the oscillatory creation and annihilation 
operators, x is the nonhnearity strength, and A = ojq—uj 
is the detuning between the mean frequency of the driv- 
ing field and the frequency of the oscillator. For f{t) ~ 1 
this Hamiltonian describes nonlinear oscillator driven by 
a monochromatic force. 

The evolution of the system of interest is governed by 
the following master equation for the reduced density ma- 
trix in the interaction picture: 

i=1.2 ^ ^ 

(3) 

where Li = ^ {N + l)ja and L2 = \/N"fa^ are the 
Lindblad operators, 7 is a dissipation rate, and N de- 
notes the mean number of quanta of a heath bath. To 
study the pure quantum effects we focus below on the 
cases of very low reservoir temperatures which, however, 
ought to be still larger than the characteristic tempera- 
ture T Tcr — h'y/kB- 

The Hamiltonian ([T]) describes wide range of physical 
systems, including nano-mechanical oscillator, Josephson 
junction device, optical fibers, quantum dots, quantum 
scissors. Part of them have been noted in section|T| Note, 
that quantum effects in a NDO with a time-modulated 
driving force including also pulsed regime have been stud- 
ied in a series of papers [46]-|51j. 

Below, we use numerical simulation of this equation 
based on quantum state diffusion method (QSD) (see, 
Ref. [55] and, for example, applications in Refs.|46|- 
[5T] and [S3]-[SS]). For clarity, in our numerical calcu- 
lation we choose the mean number of reservoir photons 
= 0. Note, that for <C 1 the above mentioned 
restriction is valid for the majority of problems of quan- 
tum optics and, particularly, for the schemes involving 
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the nano-mechanical oscillator and Josephson junction. 
In experiments, the nonlinear oscillator based on the 
current-biased J J is cooled down to T = 20mK, which 
corresponds to iV = 0.0013, whereas, Tcr = 10^^ K, for 
7 = IMHz. 

In semiclassical approach the corresponding equation 
of motion for the dimensionless amplitude of oscillatory 
mode has the following form 

^ ^-i[A + x + 2\aW]a + if{t)n~-fa. (4) 

This equation modifies the standard Duffing equation on 
the case of the NDO with time-dependent coefficient. 

III. BISTABILITY AT LEVEL OF A FEW 
QUANTA 
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FIG. 1. (Color online) The mean excitation numbers of anhar- 
monic oscillator: (1) quantum solution, (2) semiclassical so- 
lution. The parameters are as follows: A/7 = —15, x/7 = 2. 



At first, we describe the NDO driven by monochro- 
matic force [the case, f{t)—l in the Hamiltonian de- 
scribed by (1)] in bistable regime. In semiclassical ap- 
proach based on Eq. 4 the bistable dynamics is realized 
if the following inequalities are satisfied |27) : 



x(A + x) <0, 
|(A + x)/0.57| > V3, 



(A + X)' 



1.57 

A + x 



< 



7/2 
A 



X- 
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(5.) 



In this range of the parameters the typical hysteresis 
curves depending on detuning A and force amplitude 
display the system [57j. This result for the stationary 
excitation number n = |ap can be obtained by solving 
the following equation: 



(A + x + 2x|a|' 



(7/2)^ 



(6) 



and is depicted in Fig. [T] 

It is well known that whereas the semiclassical re- 
sult exhibits hysteresis-cycle behavior, the correspond- 
ing quantum mechanical result, which accounts for the 
influence of quantum noise, shows a gradual evolution. 
It is seen from the exact quantum solution for the mean 
excitation number in the following form: 
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i^(c-f l,c* + l,z) 



(A + x)2 + (7/2)2 F(c,c*,z) 



(7) 



where _F = _F2 is the generalized hypergeometric function, 
c and z are coefficients that depend on the parameters: 
c = (A + x)/x - *7/(2x) and z = 2{n/x?. 

It is also seen that the characteristic threshold behav- 
ior, determined by a drastic increase of the intensity in 
the transition region, disappears as the relative nonlin- 
earity x/7 increases. In Fig. [l] we plot both quantum 
and semiclassical solutions corresponding to Eqs. 6[7 



More detailed information on the quantum statistical 
properties of the oscillatory mode in the bistable range 
can be obtained from the analysis of the excitation num- 
ber probability distribution function p{n). In this way, 
the locations of extrema of the p(n)-function, i.e. the lo- 
cations of the most and least probable values of n, may 
be identified with the semiclassical stable and unstable 
steady states in the limit of small quantum noise level 
|57|. With the increase of x/7 the curve of locations of 
these extrema as depending on n becomes shifted from 
the corresponding semiclassical curve for the mean exci- 
tation number. 

Considering bistability at low-limit of quanta, for 
strong quantum regime, it is naturally to put the ques- 
tion what is the minimal level of excitation numbers at 
which bistable regime is displayed? Below, we discuss 
this problem on the base of quantum trajectories and 
the Wigner function to discuss quantum dynamics and 
on the Poincare section to discuss semiclassical dynam- 
ics. 

Analyzing the monochromatically driven NDO for 
strong quantum regime and in over transient time in- 
tervals, t ^ 7""'^, we use numerical method but not the 
analytical results obtained in terms of the exact solu- 
tion of the Fokker-Planck equation [57]- [59]. The rea- 
son is that the steady-state solution of the Fokker-Planck 
equation has been found using the standard approxima- 
tion method of potential equations. The validity of this 
solution has not been checked up to now in the strong 
quantum regime that requires a high nonlinearity with 
respect to dissipation. 

We use the Wigner function 



Pnm 



(8) 



in terms of the matrix elements pnm — ("IpI'ti) of the 
density matrix operator in the Fock state representation. 
Here (r, 9) are the polar coordinates in the complex phase 
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FIG. 2. (Color online) The mean excitation numbers of os- 
cillatory mode (a), the Wigner function (b), semiclassical 
Poincare section (c) , the time evolution of the excitation 
numbers along a single trajectory (d). The parameters are as 
follows: A/7 = -8, xh = 2, and n/7 = 2.7 

space plane, x = rcos0, y = rsinff, while the coefficients 
Wmnif,0) are the Fourier transform of matrix elements 
of the Wigner characteristic function. 

The properties of bistable dynamics at the level of 
a few excitation number are demonstrated on Fig. [2] 
Fig.[2]ja) shows that the mean excitation number is small 
which means that the system is operated in a deep quan- 
tum regime. Analyzing one single quantum stochastic 
trajectory for excitation numbers on Fig. [2jd) we set the 
system initially to the vacuum oscillatory state and con- 
sider time-dependence over a long time, compared to the 
characteristic dissipative time. As expected, the analy- 
sis of the time-dependent stochastic trajectories for an 
expectation number shows that the system spends most 
of its time close to one of the semiclassical bistable so- 
lutions with quantum interstate transitions, occurring at 
random intervals. 

In order to demonstrate the bistability in phase-space 
we tune the nonlinear oscillator parameters to satisfy in- 
equalities Eq. ([5| . As calculations show the Wigner func- 
tion displays two peaks (Fig. ^h)) that displays bista- 
bility, whereas Poincare section obtained from the semi- 
classical calculations of Eq. (|4| shows a single point in 
phase space (Fig. ^c)) that corresponds to the regular 
dynamics. 

Further investigation of this model allows to establish 
other properties of dissipative bistable dynamics. For 
this goal on the Fig. |3]we show bistability in phase-space 
depending on the amplitude of the external force. As we 
see for the given parameters of detuning and nonlinear- 
ity there is an intermediate range of amplitudes where 
bistability takes place. 




FIG. 3. (Color online) The Wigner functions for oscillatory 
mode in dependence on the amplitude. The parameters are 
as follows: A/7 = -8, x/l = 2, and (a) Q/j = 2.1, (b) 
Q/j = 2.3, (c) J7/7 = 2.5, (d) n/7 = 2.7, (e) fl/j = 2.9, (f) 
Q/j = 3.1. 



It is also interesting to consider behavior of the NDO 
by using the scaling properties of Eq. Indeed, it 

is easy to verify that this equation is invariant with re- 
spect to the following scaling transformation of the com- 
plex amplitude: a Aa,x xl^^-^ ^ Ail, A ^• 
^ + x(l ~ Thus, for A > 1 oscillatory excita- 

tion numbers are increased during scaling transforma- 
tion. It is interesting to analyze such scaling from the 
point of view of quantum-statistical theory and its rele- 
vance for decoherence and dissipation. Thus, using the 
scaling properties we consider the system for various ex- 
citation numbers. The Wigner functions for the scaled 
A = 2 and A = 3 are presented in Fig. |4j As we see, in- 
creasing A leads to the suppression of one of the peaks in 
the Wigner function and on the whole the bistability of 
the system is vanished. Thus, we show that such param- 
eter scaling does not occur for strong quantum regime of 
NDO. 

There is a boundary to detect bistable states at limit 
of small excitation numbers in phase-space in accordance 
with Planck uncertainty principle. Taking into account 
that AAAP > 1/2, where X and P are dimension- 
less position and momentum operators and n = |ap, 
a — X + iP, it seems that for small enough level of ex- 
citation numbers we cannot distinguish two branches of 
bistability. In this case sizes of contour plots of Wigner 
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FIG. 4. (Color online) (a) The Wigner function for A = 2, 
(b) the Wigner function for A = 3. The parameters are as 
follows: A/7 = -8, xh = 2, and n/7 = 2.7 











A- 

















(a) 



0.4 



0.2 



0.0 00 



4 




Y 



FIG. 5. (Color online) The Wigner functions showing quan- 
tum interference fringes. The parameters are as follows: 
A/7 = -8, xh = 2, n/7 = 2.7, (a) T = 0.67"^ r = 27-1, 
(b) T = 0.l7"\ and r = 27"^ 



IV. QUANTUM INTERFERENCE ASSISTED 
BY BISTABILITY 



In this section we demonstrate that it is possible to 
create quantum superposition in bistable dynamics of the 
NDO under pulsed excitation. Really, the application of 
time-dependent force can lead to transition between two 
branches of the system dynamics in the bistable regime 
and open an opportunity to generate interference be- 
tween them. However, quantum interference takes place 
for very short time intervals and disappears due to dis- 
sipation and decoherence. In order to recover the quan- 
tum interference for over transient regimes we suggest 
to apply the specific train of Gaussian pulses accord- 
ing to the model considered in section II. The results 
depicted on Fig. [5] shows that with applied Gaussian 
pulse the Wigner function has negative ranges. The 
Wigner functions describe two humps corresponding to 
bistable branches and interference pattern between them 
in phase-space. 

In order to be sure that effects of dissipation and 
decoherence in oscillatory mode for transient time are 
suppressed if the train of pulses is applied we calculate 
the purity of the state, i.e. Tr{(?). For a pure state 
Tr[p^) = 1. The results for Tr{p^) are depicted in Fig. [6] 
for over-transient regime, i.e. for time intervals, i 3> 7~ , 
in dependence on the pulse duration T and time-intervals 
between pulses r. The Fig. corresponds to the case 
of fixed T , i.e. r = 2.5. Considering the dependence from 
the duration of pulses for this case we conclude that the 
purity is maximal for very short pulses and lost with in- 
creasing of T. The opposite behavior for the fixed interval 
between pulses is realized in Fig. |6jb). 



functions are sufficiently squeezed such that two bistable 
branches are too close to each other to be distinguishable. 
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FIG. 6. The purity depending from T for fixed r — 2.67"^ (a), 
the purity depending from r for fixed T — 0.57~^ (b). The 
other parameters are: A/7 = —8, x/l ~ 2, and $1/7 = 2.7. 



V. QUANTUM DISSIPATIVE CHAOS AT 
LOW-LEVEL OF EXCITATION NUMBERS 

In this section we demonstrate that dissipative chaos is 
realized in strong quantum regime NDO on low-level of 
oscillatory excitation numbers. Chaotic regime appears 
in the NDO driven by train of Gaussian pulses, and it 
depends on the duration T of the pulses and the time 
intervals r between them. 

Many criteria have been suggested to define chaos in 
quantum systems, varying in their emphasis and domain 
of application. Nevertheless, as yet, there is no univer- 
sally accepted definition of quantum chaos. Our analysis 
is given in the framework of semiclassical and quantum 
distributions by using a correspondence between contour 
plots of the Wigner function and the Poincare section. 
Such analysis has been proposed and realized f48| in 
mesoscopic regimes of NDO with time-dependent coef- 
ficients. 

It is well known that the Poincare section has the 
form of strange attractor in phase space for dissipative 
chaotic systems while it has the form of close contours 
with separatrices for Hamiltonian systems. Thus, in this 
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FIG. 7. (Color online) (a)The Poincare section, (b)tlie contour 
plot of the Wigner function, (c)tlie excitation number time 
series with snapshots of Gaussian pulses as inset, and (d) the 
Wigner function. The parameters are as follows: xH ~ ^-^ ^ 
= 20.4, A/7 = -15, T = 0.257"\ r = 27r/57. Note that 
these distributions occur at time 'yt = 100 with the mean 
excitation number of 1.54. The range of excitation numbers 
are 1.26 to 4.98. 
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FIG. 8. (Color online) (a) The Poincare section, (b) the 
contour plot of Wigner function, (c) the excitation number 
time series with snapshots of Gaussian pulses as inset, and 
(d)the Wigner function. The parameters are: x/7 = 0-7, 
n/7 = 20.4, A/7 = -15, T = 0.2057"\ r = 27r/57. Note 
that these distributions occur at time 'yt — 100 with the mean 
excitation number of 1.74. The range of excitation numbers 
are 1.29 to 5.11. 



paper we demonstrate the quantum chaotic regime by 
means of a comparison between the contour plots of the 
Wigner functions and the strange attractors on the clas- 
sical Poincare section. In this way, we calculate the 
Wigner function in phase space by averaging an ensem- 
ble of quantum trajectories for definite time intervals. 
On the other hand, the Poincare section is calculated 
through the semiclassical distribution based on Eq. |4] 
but for a large number of time intervals: it is constructed 
by fixing points in phase space at a sequence of periodic 
intervals. Note, that such analysis seems to be rather 
qualitative than quantitative for strong quantum regime 
and the ranges of low-level oscillatory excitation num- 
bers, where the validity of semiclassical equation is ques- 
tionable. Indeed, it is shown below that semiclassical and 
quantum treatment of quantum dissipative chaos are car- 
dinally different in the deep quantum regime. 

The typical results of calculations are depicted below. 
Note, that the ensemble-averaged mean oscillatory exci- 
tation number and the Wigner functions are nonstation- 
ary and exhibit a periodic time dependent behavior, i.e. 
they repeat the periodicity of the driving pulses at the 
over transient regime. In this nonstationary regime the 
Poincare section depends on the initial time-interval tg. 
We choose various initial time to in order to ensure that 
they match to the corresponding time-intervals of the 
Wigner function. In Figs. [7j |8j [9] and [TOj we show the 
typical semiclassical and quantum distributions for the 
parameters A/7, x/7i corresponding to the chaotic 
regimes and for various durations of the Gaussian pulses. 

As we see, the figures of Poincare sections clearly in- 
dicate a classical strange attractor with fractal structure 
that is typical of chaotic dynamics. The Wigner func- 




FIG. 9. (Color online) (a) The Poincare section, (b) the 
contour plot of Wigner function, (c) the excitation number 
time series with snapshots of Gaussian pulses as inset, and 
(d) the Wigner function. The parameters are: x/7 = 0-''': 
Q/-y = 20.4, A/7 = -15, T = 0.157"\ r = 27r/57. Note that 
these distributions occur at time 'yt = 100 with the mean ex- 
citation number of 2.04. The range of excitation numbers are 
1.43 to 5.13. 

tions have spiral (helical) structures (Figs. [7|to[l0|) that 
reflect chaotic regime in analogy with the corresponding 
Poincare sections, and their contour plots are concen- 
trated approximately around the attractor. Nevertheless, 
in this deep quantum regime the different branches of the 
attractors are hardly resolved in the Wigner functions. 

It should be specified that the Wigner function for the 
regime presented in Fig. [TOjhas region of negative values. 
Obviously, this fact reflects on quantum effects in the 
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FIG. 10. (Color online) (a) The Poincare section, (b) the 
contour plot of Wigner function, (c) the excitation number 
time series with snapshots of Gaussian pulses as inset, and 
(d) the Wigner function. The parameters are: xH ~ 0-7, 
n/7 = 20.4, A/7 = -15, T = 0.l7"\ r = 27r/57. Note that 
these distributions occur at time 7f = 100 with the mean ex- 
citation number of 2.46. The range of excitation numbers are 
1.62 to 4.83. 




FIG. 12. (Color online) (a) The Poincare section, (b) the exci- 
tation number time series, (c) the contour plot of Wigner func- 
tion, and (d) the Wigner's function at the moment 7t = 100.4 
when the excitation number is minimum. The parameters are 
as follows: xh = 0.7, n/7 = 20.4, A/7 = -15, T = 0.l7"\ 
T = 27r/57. 



at this moment of the maximum excitation r 



> 




FIG. 11. (Color online) (a) The Poincare section, (b) the exci- 
tation number time series, (c) the contour plot of Wigner func- 
tion, and (d) the Wigner's function at the moment 7f = 100.6 
when the excitation number is maximum. The parameters 
are: xh = 0.7, f7/7 = 20.4, A/7 = -15, T = 0.l7-\ 
r = 27r/57. 



chaotic regime. 

The chaotic dynamics of the oscillatory mode strongly 
depends on the time-interval t. To demonstrate this point 
in Figs. 



11 and 12 we depict the Wigner function and 



the Poincare section for the oscillatory parameters used 
in Fig. [TOj however, for the other time-intervals of t 
within the duration of pulses which correspond to the 
maximal and minimal values of the number of excitation 
number. 

In the end of this section, we explore the transition 



from regular to chaotic regimes that are realized through 
varying the strength of the pulse trains by consider- 
ing cases of small excitation number. For this goal, 
we present at first the results in semiclassical approx- 
imation based on an analysis of the Lyapunov expo- 
nents of the semi-classical time series [60) . This quan- 
tity is determined L — ttt In j^ere x = 
(i?e(a), /m(a), /?), where (3 is the time variable defined 
through dl3/dt = 1 which augments Eq. ^ to create an 
autonomous system. Note that X2 and xi represent two 
trajectories that are very close together at the initial time 
to. Furthermore, At = t — to, with t — > 00. For L > 
the system shows chaotic dynamics. L — corresponds 
to the case of conservative regular systems, and L < 
indicates that the dissipative system is regular. We ex- 
amine the exponents for time intervals corresponding to 
the minimal and maximal excitation numbers of the oscil- 
latory mode in dependence from the parameter fl/j. The 
results are depicted in Figs. [13] and [14] for the constant 
parameters: A/7 and x/7- We observe a transition from 
regular to chaotic behavior at $7/7 — 12.55 for minimal 
and maximal excitation number of the oscillatory mode. 
Note that this transition occurs at the strong quantum 
regime as the excitation number ranges from the mini- 
mum n = 0.94 to the maximum n = 2.70. Thus, when 
the strength of the pulse trains fl/j is low, we observe 
regular behaviour; and as a critical threshold is crossed, 
the system behaves chaotically. Interestingly, a closer 
scrutiny of the dynamics of the system reveals a regime of 
transient chaos in the range 17.61 < O/7 < 19.56, where- 
upon the semi-classical dynamics rattle about chaotically 
for some time before settling down to regular behaviour 
which leads to a window of negative Lyapunov exponents. 
Then, beyond fl/j = 19.56, chaotic attractors are found 
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FIG. 13. (Color online) The largest Lyapunov exponent of 
the semi-classical dynamics versus the strength of the pulse 
trains. The solid (blue) curve corresponds to the moment of 
maximum excitation number at 'yt = 39.1 for 1 < Q./'y < 19 
and at yt = 39 for 19.5 < f2/7 < 26. The dashed (red) curve 
corresponds to the moment of minimum excitation number at 
"ft = 40.2 for 1 < n/-y < 8.5 and at -yt = 40.1 for 9 < Qf-y < 
26. The rest of the parameters are: x/l = 0-7, A/7 — —15, 
T = 0.l7"\ r = 27r/57. 



low-excitation, quantum noise and quantum effects play 
essential role in forming of the oscillatory dynamics, par- 
ticularly, in realization of chaotic dynamics and scenarios 
of transition from chaotic to regular regime. This state- 
ment is also confirmed by the calculation of the Wigner 
function. We observe that below the transition thresh- 
old of ri/7 = 12.55, the Wigner function has the form 
of a single hump. Beyond the transition threshold, the 
hump starts to spread and spiral begins to form. These 
changes are found to happen continuously and smoothly 
in the Wigner function as fl/j increases. However, as 
n/7 goes above 19.56, the Wigner functions are observed 
to quickly take the appearance of a strange attractor. 
Our results here thus clearly show the presence of good 
quantum-classical correspondence in the quantum and 
classical dynamical behaviour. In addition, our results 
also demonstrate that a quantum NDO, which is a form 
of quantum anharmonic oscillators, possess a rich set of 
dynamical behaviour and properties that may be poten- 
tially useful for many practical purposes [5T1 [52] . 



VI. SUMMARY 
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FIG. 14. (Color online) The plot of the maximum and mini- 
mum excitation number versus the strength of the pulse trains 
that correspond to the case of Fig. |13| The solid (blue) curve 
is for the maximum excitation number and the dashed (red) 
curve for the minimum excitation number. Note that the reg- 
ular, transient chaos, and chaos regimes of the corresponding 
semi-classical dynamics are also indicated in the plot. 



to emerge again in the Poincare section. It is important 
to note that analogous dynamical behavior are observed 
to arise at both the moments when the excitation number 
is a minimum and a maximum (see Fig. 



13) 



In summary, we have uncovered the parameters for 
which decreasing of the excitation number leads to a 
transition of the system from chaotic to regular regime 
in the classical treatment. Another situation is realized 
in the quantum treatment. Indeed, for these regimes of 



To summarize, we have studied the problem of quan- 
tum chaos and bistability at level of few excitation num- 
bers of NDO that is interacting with an external field and 
a background environment leading to dissipation and de- 
coherence. We have analyzed the formation of bistable 
behavior and chaotic regime of NDO in strong quan- 
tum regime when the ratio x/7 is chose from 0.7 to 2 
. We use a systematic numerical analysis based on nu- 
merical simulation of master equation by using quantum 
state diffusion method of quantum trajectories where we 
have varied a number of relevant oscillatory parameters: 
17/7, A/7 as well as the parameters of Gaussian pulses. 
Thus, the combination of the results provides a rather 
thorough understanding of the NDO in strong quantum 
regime as well as have presented a simple picture to un- 
derstand the formation of bistability and chaos at low- 
level of quanta in phase-space. We have also found un- 
expected features of NDO in phase-space. It has been 
demonstrated that the Wigner functions of oscillatory 
mode in both bistable and chaotic regimes realized due 
to interaction with train of Gaussian pulses acquire neg- 
ative values and interference patterns in parts of phase- 
space. We have demonstrated that in the case of bistable 
dynamics the Wigner functions describe two humps cor- 
responding to bistable branches and interference pattern 
between them in phase-space while for the chaotic regime 
the Wigner functions have spiral structures, (which cor- 
respond to strange attractor in Poincare sections), with 
deep well showing negativity of the Wigner function (Fig. 
[lO](d)). Quantum interference in phase-space is realized 
in over transient regime due to driving the oscillator by 
series of short pulses with proper parameters for effective 
reducing of dissipative and decoherence effects. Our re- 
sults can be tested with available experimental systems 
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